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Abstract 

We consider QCD with valence and sea quarks obeying different boundary conditions. We point 
out that the energy of low lying two hadron states do not depend on the boundary condition of the 
sea quarks (up to exponentially small corrections). Thus, the advantages in using twisted boundary 
conditions on the lattice QCD extraction of nucleon-nucleon phase shifts can be gained without 
the need of new gauge configurations, even in fully unquenched calculations. 
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Particle scattering is defined at infinite volumes. Hadron-hadron interactions in lattice 
QCD can be studied, paradoxically, only at finite volumes. This is because euclidean corre- 
lators with two hadrons are dominated, at large euclidean times, by the state formed by the 
two hadrons at rest. The exponential decay in imaginary time is set by the sum of the two 
hadron masses and no scattering information can be deduced from the correlator. This ob- 
servation was formalized in [1]. At finite volumes though, the hadrons are forced to interact 
and the energy levels depend on the their interaction. Thus, measuring the energy levels, 
one can deduce information about their interaction [2-4]. An elegant formula relating the 
energy levels to the phase shifts was given in [3] and we will call this approach the "Luscher 
method". 

One inconvenience with the Luscher method is that the phase shifts are obtained only at 
the discrete values of the eigenstates energies. For small boxes these values are very separated 
from each other. One can, of course, change the volume of the box and learn about the phase 
shifts at other energy values but this is, frequently, prohibitively expensive. In applying this 
method to the two-nucleon system a related problem arises. The unnaturally large value 
of their s-wave scattering length, that is, the strength of their interaction, shifts the energy 
levels far away from their non-interacting values. For box sizes smaller than about 8 fm the 
energy levels are not in the region described by effective range theory and appear at negative 
values of the energy which, from the point of view of scattering, are unphysical [5]. The 
large lattice sizes needed to extract effective range parameters values make the prospect of 
studying nuclear interaction through lattice QCD even more distant. 

In [6] it was suggested that the use of twisted boundary conditions for the quarks or, 
which is the same, simulations done under the presence of a constant background magnetic 
potential coupled to baryon number, shifts the energy levels in a calculable manner (see 
also [7]). This allows for an extra handle on the values at which the phase shifts can be 
determined. This method, however, requires the generation of new gauge configurations at 
each different value of the background field used. The purpose of this report is to point 
out that, under some circumstances, the same effect can be obtained by coupling only the 
valence quarks to the background field, which obviates the need of new gauge configurations 
generated with a background field. We also discuss how and when this method can be used 
in meson-baryon and meson-meson systems. As this paper was being finished, a similar point 
was made independently by Sachrajda and Villadoro [8], in the context of pion physics. 
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Let us consider QCD with the sea quarks satisfy standard periodic boundary conditions 
while the valence quarks have "twisted" conditions: 

q v (x,y,L) = e l0ara q v (x,y,O), 

q s (x,y,L) = q s (x,y,0), (1) 

where q v are the valence quarks, q s the sea quarks and r a = (l,r) acts on isospin space. 
This is the theory of partially twisted QCD (ptQCD) studied independently in [8]. A field 
theoretical description of ptQCD can be given using the trick used in partially quenched 
QCD [9]. We write the partition functions as 

Z = [ Dq v Dq v Dq s Dq s DqI?qe- Id ^ F2+ ^ +m ^^ (2) 

where q are "ghost" (bosonic) quarks satisfying the same boundary conditions as the valence 
quarks. Integration over the valence and ghost quark fields gives two determinants that 
cancel each other and only sea quarks appear on internal quark loops. The ghosts violate 
the spin-statistics theorem and ptQCD is not a unitary theory. The "wrong" statistics for 
the ghosts will generate a Hilbert space with a non-positive metric. Notice that the isospin 
(and baryon number) of valence, sea and ghost quarks is separately conserved. 

We can trade the twisted boundary conditions by the presence of a constant gauge field 
potential by performing a field redefinition similar to a gauge transformation, except for 
being discontinuous at z = L: 

■ zB a 

q v (x,y,z) -> e l ~ ra q v (x,y,z), 

■ zO a 

q(x,y,z) -> e l — ra q(x,y,z). (3) 

(4) 

The effect of this transformation is to change the lagrangian in Eq. (2) to 

C q = q v (fl + iA + m)q v + q(p + ifi + m)q + q s {Ip + m)q s , (5) 

with = (0, 6 a T a /Lz). Obviously, this lagrangian has now a SU isospin (2) x £7^(1) gauge 
symmetry of the form 
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q v -> U(x)q v , U(x) G SU(2) x C/(l) 

g„ -> U(x)q v 

q s -> g s 

A„ -> [/(^AC/^x) -iUd^ix). (6) 

The maximum value of the external field A is A = tt/L which, in all interesting cases, is 
of the order of m n or less. Larger values of A can be eliminated by a (continuous) gauge 
transformation and has no physical consequence. 

The sea sector of the theory is identical to QCD since the valence and ghost determinants 
cancel and there is no back-reaction of the dynamics of the valence quarks on dynamics of 
the sea quarks. That means that the condensate (q s q s ) is the same as in QCD. Due to 
the symmetries of the theory (q s q s ) = (q v q v ), as discussed by [10]. The whole pattern of 
symmetry breaking is described by the diagram below: 

SU L (2N\N) x SU R (2N\N) x U l+r {1) > SU L+R (2N\N) x U L+R (1) 

m^O 

SU L+R (™\N) x U L+R {1) 

SU L+R (N\N) x SU L+R (N) x U{1) x U(l) { ===^SU L+R (N\N) x SU L+R (N) x U(l) x U(l) 

In the upper left corner we have the full symmetry in the absence of either quark masses 
or external field. On the lower right corner the remaining symmetry after spontaneous 
chiral symmetry breaking, a common quark mass for all quarks and the external field are 
included. Double arrows indicate spontaneous breaking, single arrows the explicit breaking 
of symmetries. 

This pattern of symmetry breaking implies the existence of Goldstone bosons made up 
of valence, sea and ghost quarks (in addition to the ones formed only by valence quarks) 
forming multiplets of the remaining symmetry SUl+r(N\N) x SUl+r(N) x {7(1) x {7(1). In 
the N = 2 case we have not only pions made of u and d quarks but also ghost pions made of u 
and d ghosts, sea pions made of u s and d s , and mixed combinations including fermionic states 
like, for instance, u'jsd, ■ ■ ■ . A similar phenomena occurs with other hadrons. In hadronic 
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low energy effective theories, hadrons containing sea and ghost quarks occur inside loops 
and serve the purpose of canceling the contribution of valence quark loops and substituting 
them by the correct ones with sea quarks instead. The construction described above is 
very similar to the one used in partially quenched QCD for the meson [10-13], one baryon 
[14, 15] and two-baryon [16] sectors. An alternative (and equivalent) formalism is given by 
the replica method [17, 18]. 



A. Nucleon-nucleon scattering 

Since the number (and isospin) of valence quarks is conserved, all intermediate states in 
two-nucleon processes contain the two initial baryons. For the sake of argument we first 
consider the regime Q <C m n , where Q is the typical momentum of the state considered and 
the pion mass. The effective theory in this regime is described by 



C NN = N\iD + ^)A - Co (A T PA) f (N T P t N) + ^C 2 [(A T P,A)t (N T O t N) + h.c] +■■■ 
+ N*(iD + ^)A - C (N T P t N)^N T P t N + ■■■ 

+••• (7) 

where the dots represent higher derivative terms, A is the valence nucleon field, A the 
field of a nucleon containing one sea quark and Pi is a spin-isospin projector and Oi = 
PjD 2 + D 2 Pj — 2DPjD. Also not shown explicitly are the terms containing the nucleons 
with other combinations of valence, sea and ghost quarks. The background field can only 
appear within covariant derivatives DA = DN + iAN since the underlying theory has 
a gauge symmetry (Eq. 6). There are two important points to notice. The first is that 
the conservation of valence isospin and quark number forbids terms coupling valence and 
non- valence quarks so the valence nucleons effectively decouple from its analogues containing 
ghost and/or sea quarks. The second one is that the coefficients in Eq. 7 encapsulate physics 
in the scale ~ l/m n and have only an exponentially small dependence on the volume and 
boundary conditions of order e~ m7rL . 

We can obtain the Luscher formula by comparing two calculations in this effective theory, 
the scattering amplitude in infinite volume and a finite volume two-nucleon correlator. They 
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FIG. 1: Graphs contributing to nucleon-nucleon scattering in the Q <C m n regime. 

are both given by the iteration of the same short distance interaction, as shown on Fig. 1. 
The correlator can be computed exactly (using dimensional regularization), to all orders in 
the low-energy expansion, as shown below: 



C(E,p) = J dte iEt (0\TN T (p + p o ,t)PN(-p + p o ,t) . N^-p + p , 0)P f iV*(p + p , 0)|0) 

~ E E -E n + iO l{En |iVt( " P + P °' °) ptiV *(P + P°' °)l°>| 2 
q 

1 
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where the sum is over the allowed momenta 

Qx,y — -^ n x,y) 

g 2 = T (n 2 + ^), (9) 

with n x ,n y and n z integers, and po = <f>/Lz is the minimum momentum allowed in the 
lattice. 

By the other hand the same combination of low energy constants appearing in Eq.(8) 
also determines the scattering in the infinite volume limit 



l-/oE„^(/i)A: 2 « 
4tt 1 

~ M k cot 5 - ik ' ^ ' 

Comparing Eqs. 8 and 10 we can relate the shift in the pole of the correlator in Eq.(8) to 
the infinite volume phase shifts. This provides the generalization of the Luscher formula to 
the twisted boundary condition discussed in [6]. 

This shows that in ptQCD, baryon-baryon scattering at very low energies is the same as 
in QCD with twisted (valence and sea) quarks. The numerical estimates of the influence of 
the boundary conditions on the two-nucleon states is then unchanged from [6]. 
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This result is valid up to a higher energy regime. In fact, it is correct up to the threshold 
for pion production. To show this we need to modify the derivation given above. At higher 
momenta of order Q ~ m n the appropriate effective theory contains mesons explicitly. 
Valence baryons can then couple to non- valence baryons through the exchange of non- valence 
mesons. We should distinguish two kinds of scattering graphs: irreducible ones (which do 
not come apart by cutting two valence nucleon lines) and reducible ones. The non-valence 
nucleons and mesons appear inside the irreducible parts of diagrams only (see e.g. Fig. 2), 
the finite volume effects of which are of order e ~ mirL . It remains true that intermediate 
states, sensitive to finite volume effects, are always composed of two valence nucleons. The 
full amplitude results from iterating the irreducible graphs or, what is the same, solving the 
Lippman-Schwinger equation: 

T(k,p)=V(k,p)-J ^T(k,q) E _£ +iQ V(q, P ), (11) 

where V(k,p) is the sum of the irreducible graphs (potential). The scattering amplitude is 
related to the /C-matrix satisfying 



lC(k,p)=V(k,p)- 




through 

r(M = (i + ^m*)) ■ (is) 

The /C-matrix is related to the phase shifts by l/]C(k, k) = Mk cot 5 /An, as can be seen from 
Eqs. 10 and 13. Eq. 12 is modified at finite volume only by the change of the integral over 
intermediate states by a sum 

f d 3 q 1 v-^ . . 27T 2tt 6 

y ^•••-^^•••,with q x , y = -n x , y ,q z = -(n z + -). (14) 

q 

The typical momentum of this sum is determined by the scale of V(k,p), in our case, 
q ~ m % . As long as the box satisfies m n L 3> 1 the effect of the discretization of the 
intermediate momenta is small and the /C-matrix at finite volume differs from its infinite 
volume counterpart by terms proportional to e~ m7vL . For instance, the leading source of 
boundary condition dependence in nucleon-nucleon scattering is the one-pion exchange part 
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of the nuclear potential, as this is the longest range contribution of the potential. The 
equation determining the /C-matrix at finite volume is (after projecting in the spin singlet) 



E(M = nk,p)-^IC(k,y) ^ (15) 

q ^ M yH ^' 

= V(M-^E / d\6(i- v)Hk, y) E \ ^ {y _ fy + m , (16) 



mI ^ <1 Ae-'™- L (4+2cos(</0) 



(17) 



where the allowed values of q are like in Eq. (9) and I = 2irn/L. The factor 4 + 2cos(0) 
comes from summing e~ tlz ^ over the six points in wave vector space with unit length. Notice 
that this is not true for the T- matrix. In the infinite volume limit the integral in Eq. 11 is 



dominated by q ~ V ME and for these values the integral is not well approximated by the 
discrete sum. The argument presented above in the Q <C m n regime can now be applied to 
the ^/Mm^ > Q ~ m T regime just by changing J2 n C 2n (ME) n by 1/JC(VME, y/ME). 
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FIG. 2: Contribution to nucleon-nucleon scattering involving sea quarks. 



1. Meson-meson and meson-nucleon scattering 

In systems containing mesons, the presence of valence aniz-quarks may invalidate the 
point in the previous section. It is important to distinguish two kinds of channels. In the 
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first we have those channels where the initial valence quarks cannot annihilate with each 
other, for instance, the two-pion 1 = 2 channel or the 1 = 3/2 pion-nucleon channel. In 
this case, intermediate states have to include the original mesons and sea quarks can appear 
only on t-channel exchanges (see Fig. (3)). 




FIG. 3: Diagrams contributing to I = 2 ir — ir scattering in QCD (top row), in an effective theory 
(%PT) valid below the confinement scale (second row) and in a effective theory valid below the 
pion mass (lower row). Sea quarks and mesons containing sea quarks are denoted by red, dashed 
lines. 

The situation in these channels is identical to the two-nucleon case discussed above. The 
non-unitarity due to the different boundary conditions for sea and valence quarks is relatively 
harmless at low energies and manifests itself only on terms suppressed by e~ mnL . The end 
result is that the results of ptQCD are the same as the ones in fully twisted QCD, and one 
can relate their energy levels to QCD phase shifts through the generalization of the Luscher 
formula discussed in [6] . In the two-pion, 1 = 2 channel, the external field needs to couple to 
the third isospin component (or alternatively, we can use opposite boundary conditions for 
up and down quarks), if it is to have any effect on the energy levels. This coupling breaks 
isospin but conserves /3, and preventing mixing between the I = 2, Is = 2 state from the 
1 = state. 

In systems where the incoming valence quarks can annihilate the situation is more com- 
plicated. Intermediate states with only hadrons containing at least one sea/ghost quarks 
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FIG. 4: Diagrams contributing to / = ir — tt scattering in QCD (top row), in an effective theory 
(%PT) valid below the confinement scale (second row) and in a effective theory valid below the 
pion mass (lower row). Sea quarks and mesons containing sea quarks are denoted by red, dashed 
lines. 

can occur (see Fig. (4)). In these diagrams the rest mass of the incoming hadrons is released 
and the intermediate states can go on-shell. They generate cuts in the amplitude all the 
way down to zero momentum and cannot be integrated out. Even at low energies the lack 
of unitarity is evident and there is no way to relate energy levels to phase shifts in an exact 
way. One can, of course, develop a chiral perturbation theory adequate to partially twisted 
QCD in the molds of partially quenched chiral perturbation theory, and then relate lattice 
observables to low energy constants . This method, however, can only be accurate to a 
certain order in the chiral expansion and it is unclear if there is any advantage in using 
partially twisted QCD in these cases. 
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Science Council of ROC. 
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